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We show that it is possible to reconcile NMR and neutron scattering experiments on both 
La2-xSrxCu04 and YBa2Cu306+x, by making use of the Millis-Monien-Pines mean field phe- 
nomenological expression for the dynamic spin-spin response function, and reexamining the standard 
Shastry-Mila-Rice hyperfine Hamiltonian for NMR experiments. The recent neutron scattering re- 
sults of Aeppli et al |l| on Lai.86Sro.i4Cu04 are shown to agree quantitatively with the NMR 
measurements of ^^T\ and the magnetic scaling behavior proposed by Barzykin and Pines, jij The 
reconciliation of the ^^O relaxation rates with the degree of incommensuration in the spin fluctuation 
spectrum seen in neutron experiments is achieved by introducing a new transferred hyperfine cou- 
pling C' between ^^O nuclei and their next nearest neig hbor Cm2+ spins; this leads to a near-perfect 
cancellation of the influence of the incommensurate spin fluctuation peaks on the ^^O relaxation 
rates of La2-xSrxCu04. The inclusion of the new C' term also leads to a natural explanation, 
within the one-component model, the different temperature dependence of the anisotropic ^^O re- 
laxation rates for different field orientations, recently observed by Martindale et al. The measured 
significant decrease with doping of the anisotropy ratio, ^^R —^^ Tiab/''^Tic in La2-xSrxCu04 sys- 
tem, from ^^R = 3.9 for La2Cu04 to R ~ 3.0 for Lai.85Sro.i5Cu04 is made compatible with the 
doping dependence of the shift in the incommensurate spin fluctuation peaks measured in neutron 
experiments, by suitable choices of the direct and transferred hyperfine coupling constants Ap and 
B. 



I. INTRODUCTION 



The magnetic behavior of the planar excitations in the cuprate superconductors continues to be of central concern 
to the high temperature superconductivity community. Not only does it provide significant constraints on candidate 
theoretical descriptions of their anomalous normal state behavior, but it may also hold the key to the physical origin 
of high temperature superconductivity. Recently two of us have used the results of NMR experiments to determine 
the magnetic phase diagram for the La2-xSrxCuOi and Y Ba2Cu30Q+x systems. [|[ We found that for both systems 
bulk properties, such as the spin susceptibility, and probes in the vicinity of the commensurate antiferromagnetic 
wave vector (7r,7r), such as ^^Ti, the ^^Cu spin relaxation time, and ^^T2g, the spin-echo decay time, display z — 1 
scaling and spin pseudogap behavior over a wide regime of temperatures. On the other hand, the neutron scattering 
experimental results of Aeppli et al on Lai^^^Sro iA^CuOn which probe directly x"((7,w), the imaginary part of 
the spin-spin response function, while supporting this proposed scaling behavior, at first sight appear incapable of 
explaining NMR experiments on this system. 

This apparent contradiction between the results of NMR and neutron scattering experiments, both of which probe 
x(q, w) in Lai.^^SrQuCuOA, is but one of a series of such apparent contradictions. For example, in the Y Ba2Cu-iOQ+x 
system, NMR experiments on ^^Cu and ^^O nuclei in both Y Ba2CuzO'j [|| and Y Ba2Cu30QS'i (|] require the presence 
of strong antiferromagnetic correlations between the planar Cu^^ spins, and a simple mean field description of the 
spin-spin response function with a temperature dependent magnetic correlation length ^ ^ 2, was shown to provide 
a quantitative description of the measured results for ^^T and ^^Ti in Y Ba2Cu^O-!, Q and Y Ba2Cu^OQ,Q-i. 
Yet neutron scattering experiments on YBa2Cu^Oj ]8|-p^ and Y Ba2Cu-iOQ,Q?,, 10 find only comparatively broad, 
temperature-independent, peaks in x" {li corresponding to a quite short ^ 1) temperature-independent magnetic 
correlation length. The apparent contradiction is especially severe for the La2-xSrxCu04 system, where neutron 
scattering experiments show at low temperatures four incommensurate peaks in the spin fluctuation spectrum, whose 
position depends on the level of Sr doping, [ p^ while the quantitative explanation (using the same one-component 
phenomenological description which worked for the YBa2Cu306+x system) of the measurements of ^'^Ti and ^^Ti 
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in this system requires that the spin fluctuations be peaked at (tTjTt), or nearly so. Viewed from the NMR 

perspective, there are two major problems with four incommensurate spin fluctuation peaks. First, the Shastry- 
Mila-Rice (SMR) form factor, (l^Jl^ which, provided the peaks are nearly at (7r,7r), effectively screens neighboring 
^^O nuclei from the presence of the strong peaks in the nearly localized Cv?'^ spin spectrum required to explain the 
anomalous temperature-dependence behavior of ^'^Ti , fails to do so for the considerable degree of incommensuration 
in the peaks at (tt, [7r±5]) and {[k±5],ii) seen in Lai,%QSrQ,i/^CuOi. | p^ , p7 18| As a result ^''Ti picks up a substantial 



anomalous temperature dependence which is not seen experimentally. Second, with the doping-independent values of 
the hyperfine couplings which appear in the SMR form factors for a commensurate spectrum, the calculated anisotropy 
oi^^Ti for the incommensurate peaks seen by neutrons is in sharp variance with what is seen in the NMR experiments. 



Two ways out of these apparent contradictions have been proposed. One view is that the spin fluctuation peaks 
seen in the neutron scattering experiments reflect the appearance of discommensuration, not incommensuration; on 
this view, the La2_xSrxCu04 system contains domains in which the spin fluctuation peaks are commensurate (so that 
there are no problems with ^^Ti), but what neutrons, a global probe, see is the periodic array of the domain walls. 

A second view is that a one-component description of x(q, w) is not feasible; rather, the transferred hyperfine 
coupling between the nearly localized Cv?^ spins and the ^''O nuclei is presumed to be very weak, and the ^^O nuclei 



are assumed to be relaxed by a different mechanism, whence the nearly Korringa-like behavior of ^'^Ti. |18| A further 
challenge to a one-component description has come from the very recent work of Martindale et al, ^ who find that 
their results for the temperature-dependence of the planar anisotropy of ^^Tiq, for different field orientations appear 
incompatible with a one-component description. 

In the present paper we present a third view: that the one-component phenomenological description is valid, but 
what requires modification are the hyperfine couplings which appear in the SMR Hamiltonian which describes planar 
nuclei coupled to nearly localized Cu^^ spins. We find that by introducing a transferred hyperfine coupling C", 
between the next nearest neighbor Cu^^ spins and a ^^O nucleus, the nearly antiferromagnetic part of the strong 
signals emanating from the Cv?^ spins can be far more effectively screened than is possible with only a nearest neighbor 
transferred hyperfine coupling, so that the existence of four incommensurate peaks in the La2-xSrxCu04 system can 
be made compatible with the ^^Ti results. We also find that by permitting the transferred hyperfine coupling, _B, 
between a Cv?^ spin and its nearest neighbor ^^Cu nucleus to vary with doping, we can explain the trend with doping 
of the anisotropy of ^^Ti in this system. We then use these revised hyperfine couplings to reexamine the extent to 
which the recent results of Aeppli et al [Q on Lai.86Sro.i4Cu04 can be explained quantitatively by combining the 
Millis-Monien-Pincs (hereafter MMP) response function Q with the scaling arguments put forth by Barzykin and 
Pines. 1^ We find that they can, and are thus able to reconcile the neutron scattering and NMR experiments on this 
member of the La2-xSrxCu04 system. 

We present as well the results of a reexamination of the NMR and neutron results for the YBa2Cu306+x system. 
Here we begin by making the ansatz that it is the presence of incompletely resolved incommensurate peaks which is 
responsible for the broad lines seen in neutron experiments. We follow Dai et al, ^ who suggest the increased line 
width for YBa2Cu3 07 seen along the zone diagonal directions reflects the presence of four incommensurate peaks, 
located at Qi = (tt ± ^, tt ± ^ ) , a proposal which is consistent with the earlier measurements of Tranquada et al for 
YBa2Cu30e,Q [Q. We then find that incommensuration can be made compatible with NMR experimental results 
provided the transferred hyperfine coupling constant, B, is doping dependent in this system as well. Moreover, on 
considering ^''Ti for YBa2Cu307, we find that the anomalous temperature dependence of the planar anisotropy of 
^^Ti measured by Martindale et al ^ constitutes a proof of the validity of our modified one-component model. Thus 
incommensuration combined with the presence of the next nearest neighbor coupling, C", leads to results which are 
consistent with experiment, and we are able to preserve the one-component description of the planar spin excitation 
spectrum. 

The outline of our paper is as follows: In Section II we review the SMR description of coupled Cu^~^ spins and nuclei 
as well as the mean field description of x(q, w), and examine the modifications brought about by incommensuration 
and next nearest neighbor coupling between Cu^~^ spins and a ^^O nucleus. In Section III wc review the experimental 
constraints on the hyperfine coupling parameters, and present our results for their variation with doping in both the 
La2-xSrxCu04 and YBa2Cu306+x systems. We show in Section IV how the ^^Cu NMR results can be reconciled with 
neutron scattering results on Lai.86Sro.i4Cu04, while in Section V we present a quantitative fit to the ^^Tic results 
for the iai.85S'ro.i5Cu04 based on the four incommensurate peaks in the spin fiuctuation spectrum expected from 
neutron scattering. We show in Section VI how the anomalous results of Martindale et al ^ for the YBa2Cu306+x 
system can be explained using our modified one-component model, and in Section VII we present our conclusions. 
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II. A GENERALIZED SHASTRY-MILA-RICE HAMILTONIAN 



On introducing a hyperfine coupling C'^ ^ between the the ^^O nuclei and their next nearest neighbor Cu^"^ spins, 
we can rewrite the SMR hyperfine Hamiltonian for the ^^Cu and ^^O nuclei as: 



nn 

Hh f = (r,) [ Ac...i3Si3 {r,) + BY,Sa {vj ) 



(1) 



where ^a,/3 is the tensor for the direct, on-site coupling of the ^^Cu nuclei to the Cu^+ spins, B is the strength 
of the transferred hyperfine coupling of the ^^Cu nuclear spin to the four nearest neighbor Cu^+ spins, Ca,[3 is the 
transferred hyperfine coupling of the ^^O nuclear spin to its nearest neighbor Cu^"^ spins, and ^ its coupling to 
the next nearest neighbor Cu^+ spins. The indices "nn" represent nearest neighbor electron spins to the specific 
nuclei, "nnn" the next nearest neighbor Cu^^ spins. As we shall see below, inclusion of the ^ term enhances the 
cancellation of the anomalous antiferromagnetic spin fluctuations seen by the "'^^O nucleus, and therefore reduces the 
leakage from incommensurate spin fluctuation peaks to the ^^O relaxation rates. It thus enable us to reconcile the 
measured ^^O relaxation rates with the neutron scattering experiments for both La2-xSrxCu04 and YBa2Cu306+x- 
The spin contribution to the NMR Knight shift for the various nuclei are: H 
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Kb = 



{A, + 4B)xo 
2{Cp + 2C^)xo 
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(2) 



Here, we have incorporated the new term into the ^^O Knight shift expression for ^^if/j, while the others remain 
their standard form as in Ref. 7„ are various nuclei gyromagnetic ratios, 7e is the electron gyromagnetic ratio, 
and xo the static spin susceptibility. The indices c and ah refer to the direction of the applied static magnetic field 
along the c-axis and the a6-plane. The spin-lattice relaxation rate, ("Ti)^^, for nuclei a responding to a magnetic 
field in the (3 direction, is: 



arji—l 



2^ 



q 



^/3(q)x"(q,w^o) 



(3) 



where the modified SMR form factors, "F/3(q), are now given by: 

^^Fc = [Aab + 2B{cos q^a + cos qya)f 
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Fa = 



[Ac + 2_B(cos qxa + cos %a)] 
"2 



2 E cos2^(C„, +2C;^cosg,a)', 



(4) 



Here, a' and a" are the directions perpendicular to a. The form factor ^^F^^^ is the filter for the ^^Cu spin-echo 
decay time ^^Tsg @]: 
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^2G ~ 



■69 l^Yr^" 



(q)'[x'(q,o)]^ 



^E^aY^(q)x'(q'0) 



(5) 



The values of the hyperfine constant Cq and can be determined by the various ^^O Knight shift data. In fact, 
we may obtain these new values from the "old" values of the hyperfine coupling constant, C^'^, which have been well 



3 



established from fitting the Knight shift data. Note we use Cq. to represent the new nearest neighbor hyperfine 
couphng constant, while the old hyperfine coupling constant is written explicitly as C^'^ throughout the paper. In 
order not to change the Knight shift result of the previous analysis the new hyperfine coupling constants should 
satisfy the following requirement: 

Co. + 2C'^ = = CaCf" (6) 

where Cq = C^'^/C°^'^, and c denotes the case of a magnetic field along the c-axis. For YBa2Cu306+x, from the 
previous analysis of Yoshinari and Martindale et al we have for a field parallel to the Cu-0 bond, C||=l-42, 
and — 0.91 for a field perpendicular to the Cu-0 bond direction, while Cc = 1- On introducing = C^/Cc we 
obtain: 



2rc + r 

C'.-Cf^^^ (7) 
Substituting these values of Ca and C'^ into Eq.(Q), we obtain the new ^^O form factor in terms of C"'"^: 

"^"^Tiriry^ ^ cos2^[CQ.(l + 2r,)-2r„,+2r„.cosg,a]2 (8) 

ai—a'^a" 

Although may well be anisotropic (as Cq is), in the absence of detailed quantum chemistry calculations, (which 
lie beyond the purview of the present paper) we assume to be isotropic for illustrative purposes, in which case 
r±_ — 7'|| ~ Tc ^ r = C jCc- In Fig.l, we compare our modified form factor Fc, Eq(H), with the standard 
SMR form. It is seen that with a comparatively small amount of next nearest neighbor coupling, corresponding to 
r = C jCc = 0.25, the new form factor is reduced significantly near (tt/o, 7r/a), and is some 30% narrower near q = 0. 
This indicates that the oxygen (^^TiT)~^ is less likely to pick up the anomalous antiferromagnetic contribution near 
(tt/c, 7r/a), even when the anomalous spin fluctuation is slightly spread away from (tt/o, tt/o). 

We adopt the phenomenological MMP expression for the spin-spin correlation function, modified to take into 
account the presence of four incommensurate peaks at Qi near (7r/a,7r/a) Ql, 

xq,^ =TZ^-r— 7 — TTxm — —I — + 1 — ■ — ^ 

4 ^ 1 + (q — y^iYi^ — luo/bJsF 1 — iTTw/l 

Here the first term, often called xaFi represents the anomalous contribution to the spin spectrum, brought about by 
the close approach to antiferromagnetism of the Fermi liquid in the vicinity of the peaks at q = Qi determined by 
neutron scattering experiments p],p3[. For Lai.86Sro.i4Cu04, Qi = (tt/o, [7r±(5]/a), ([tt ± (5]/a, 7r/a), with 5 = 0.2457r. 
In Eq.(|), bJSF is the characteristic frequency of the spin fluctuations, ^ is the correlation length, and a is the scale 
factor (in units of states/eV, where /is is the Bohr magneton), which relates xQ; to thus the height of each of the 
four peaks is, 

XQ. = f CVI- (10) 

The second term on the right-hand side of Eq.(^, usually called xfl-, is a parameterized form of the normal Fermi 
Liquid contribution, which is wave- vector independent over most of the Brillouin zone; F is of order the Fermi energy. 
The static bulk susceptibility xoi which is generally temperature dependent, has been determined for La2-xSrxCu04 
and YBa2Cu306+x from copper and oxygen Knight shift experiments. ||^ For a system with any appreciable antifer- 
romagnetic correlations (^ ^ a), the normal Fermi liquid contribution is small compared to xaf for wave vectors in 
the vicinity of Qi, and plays a negligible role in determining (^^TiTc)"^; however, because of the filtering action of 
-^''fo,, it makes a significant contribution to {}'^TiT)~^. Note that because the MMP expression for xaf is a good 
approximation only for wave vectors in the vicinity of the antiferromagnetic wave vector Qj, the above expression 
should not be used in calculating long wavelength properties, such as the Knight shift of ^^O. 

For the frequently encountered case of long correlation lengths (<^ > 2a), in calculating the various ^^Cu relaxation 
rates one can approximate x"(q, a;) by x"(Qi, i^)'^(q — Qi)- O^^^ can then replace Eq.(|3|) and (||) by the following 
analytic expressions: 
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(11) 

(12) 



Another important quantity, the anisotropy ratio of the ^'^Cu spin-lattice relaxation rates, which has been measured 
for La2-xSrxCu04 and YBa2Cu306+x at various doping concentrations, provides a direct constraint on the hyperfine 
coupling constants, and B. For ^ 2, this anisotropy ratio, ^^R can be written as. 
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For the case of La2_xSrxCu04, where the peaks are located at = (vr/a, [tt ± S\/a)^ ([tt ± S\/a, ir/a), we then have 



'R 



Ac - 2B(1 + cos(5) 
Aab - 2B{l + cosS) 



(14) 



For YBa2Cu30g+x; as indicated in the Introduction, on assuming the broad (vr/a, tt/o) peak seen in neutron scattering 
experiments ||^-[ll|] reflects the presence of four unresolved overlapping incommensurate peaks located along the zone 
diagonal directions, we may write 



and the anisotropy ratio becomes. 
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Q, = ([7r±J]/a, [7r± 6]/a), 
{Ac -4:B cos "I 



{Aab ~4.B cos S)\ 



(15) 



(16) 



Numerical calculations of the ^^O relaxation rates show that these rates can deviate significantly from those obtained 
by approximating the Xaf ^ Qi) function. We therefore calculate the ^^O relaxation rates numerically, using 
Eqs.(|) and (I). 



III. THE DIRECT AND TRANSFERRED HYPERFINE CONSTANTS 

Seven years of NMR experiments on aligned powders and single crystals of the cuprates have produced a signif- 
icant number of constraints which must be taken into account in selecting the hyperfine constants which enter the 
SMR Hamiltonian. Thus experiments which determine the ^^Cu nuclear resonance frequency in the AF insulators, 
YBa2C'u30e Q and La2CuOi yield similar results for the product of {AB — Aab) and ^e//, the effective moment 
of the localized Cu^+ spins, j2^ 

^lef f (4B -Aab)^ 79.65 ± 0.05 kOe (YBazCuaOe) (17) 

^J.eff{4:B - Aab) = 78.78 kOe (La2Cu04) (18) 

On using the value, fieff = 0.62/is, determined by Manousakis |2^ for the 2D spin 1/2 Heisenberg antiferromagnet, 
we then find 

AB - Aab = 128.5 kOe/^lB (YBa2Cu306) (19) 

AB - Aab = 127 kOe/fiB (La2Cu04) (20) 

A second set of constraints comes from ^^Cu Knight shift experiments. To a high degree of accuracy, in the 
YBa2Cu306-(-x system the ^^Cu Knight shift in a magnetic field along the c-axis is temperature independent in both 
the normal and superconducting state, and hence reflects only the chemical shift. The absence of a spin contribution 
means that for this system, 

Ac + AB:^0, (21) 
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independent of doping level. A third set of constraints is obtained from measurements of the anisotropy of the ^^Cu 
spin-lattice relaxation rates; for Y Ba2CuzOj one finds ^^R = 3.7 ± 0.1. To the extent that Aab, A^, and B are 
independent of doping level in YBa2Cu306+x, and the spin fluctuation peaks are commensurate (or nearly so) for 
this system, one then finds from Eqs. (|6|), ([l^), and (pT]), that 

B = 40.8 kOe/nB 
Ac = -163 kOe/fiB 

Aab^M kOe/^iB (22) 

in agreement with the analysis of Monien, Pines, and Takigawa. |^] These values are consistent with the constraint 
on (4B -|- Aab) obtained by Ishida et al for Lai.s^Sro i^CuOi, from the slope of a plot of their direct measurement of 
Xo(,T) against their measured value of ^^Kab{T), they found |p9|| 

iB + Aab = 189 kOe/^iB- (23) 

It seemed natural therefore to conclude that not only were Aab, A^ and B independent of doping for the YBa2Cu306+x 
system, but that the corresponding values for the La2_xSrxCu04 system were likewise doping independent and were 
virtually identical with those deduced for YBa2Cu3 06+x- 

If, however, the spin fluctuation peaks in the La2_xSrxCu04 system are incommensurate, the assumption that the 
hyperfine constraints for this system are doping independent is no longer tenable for this system, as may be seen by 
comparing the measured values of ^^R for the La2_xSrxCu04 system shown in Table I with the values calculated 
using Eqs(p2|), and using the doping dependence of the degree of incommensuration determined in neutron scattering 
experiments |]l^ , 6 ^ 1.75a;, where x is the Sr doping level. As may be seen in Table I, the calculated trend with 
doping is opposite to that seen experimentally. Since the quantum chemical environment responsible for the direct 
hyperfine interaction Aa is not expected to vary substantially with doping, the most likely culprit in Eqs.(p^ is the 
assumption that the transferred hyperfine coupling constant does not vary appreciably with doping; indeed, if B 
increases sufficiently rapidly with doping, with Aab and Ac fixed, one can find a doping dependence of ^"^i? which 
is more nearly in accord with experiment. This means abandoning for the La2_xSrxCu04 system the constraint, 
Ac ~ —4B, which works so well for the YBa2Cu306+x system. 

Suppose then one starts anew with the insulator, La2CuOi. On making use of Eqs. ( p^ ) and (|2|) and taking 
^^R = 3.9, in accord with the result of Imai et al |3^] at 475K, one finds readily that 

Aab-Ac = 2Q?, kOe/fiB- (24) 

On turning next to Lai.ssSro i^CuO^, taking ^"^i? — 3.0, in accord with the recent measurement of Milling and 
Slichtcr, [ pT| using the result of Ishida et al, ||29t| Eq.(^), and assuming that Aap is independent of doping, one then 
finds B — A8 kOe/iiB and Aab — ~3 kOe/fiB- This result is, however, unrealistic. A straightforward calculation 
using the expressions adapted by Monien et al from the work of Bleaney et al |32j , 

4 62 

Ac^395[-k~--—j] kOe/^iB 

Aafc = 395hK-^-y7] kOe/^lB (25) 

In Eqs.(p5|), 7 = X/E^-y is the dimensionless ratio of the spin-orbit coupling for a Cu^+ ion, A ~ — 710cr7i^^, to the 
excitation energy from the ground state of the ^'^Cu d^2_y2 orbital of the various ^^Cu d states, E^y ^ E^z ^ Eyz ^ 
2eV; with these typical values, 7 — 0.044 ± 0.009; (^) which enters as a multiplicative factor in Eq. ( p5| ) is taken to 
be 6.3a-3. With the value of 7 = 0.0469 obtained using Eq. (H), 

= (-395k -f 142) kOe/fiB- (26) 

On taking the core polarization k = 0.26 ± 0.06 we then get, for k in the vicinity of its plausible upper limit, 
0.32, 

Aab > 16 kOe/fiB- (27) 

In order to satisfy the above constraints, we next assume that the anisotropy, ^^R, for Lai.85Sro.i5Cu04 is at 
the upper end of the range quoted by Milhng and Slichter, and take ^^R = 3.2; we next take Aab = 18 kOe/fiB 
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(corresponding to k = 0.316), a value close, but not at, the estimated minimum value for Aab- We then have, from 
Eq. (p^, Ac = -185 kOe/^iB and, from Eq.jl^) for ^^R, i?o.i5 = 51 kOe/^B, while for the insulator, we find from Eq. 
(pO|), Bo — 36.1 kOe/^B- With these hyperfine constants we find for Lai.85S'ro.i5Cu04 that AB + Aab = 222 kOe/jis, 
some 17% above the value obtained by Ishida et al while for this system, the ratio of the spin contributions to 
the Knight shift for fields parallel and perpendicular to the c-axis is 

e^Kab ^B + Aab 

The slight temperature variation of ^^Kc which follows from this choice of parameters would not be detectable, 
consistent with the measurements of Ohsugi et al. |Q 

For intermediate levels of Sr doping, if we assume that the change in B induced by doping scales with the doping 
level, we obtain the results for Lai gSro ^CuOi and Lai,s6Sr()nCu04 given in Table III. Also given there are the 
corresponding results for ^^Ti and ^^T2g and related quantities of interest in analyzing NMR experiments. We note 
that to obtain ^^R = 3.5 for Lai gSro iCuO^, one needs a transferred hyperfine coupling, B = 37.8 kOe/fis, which 
is considerably lower than that obtained by direct interpolation. 

We turn next to the YBa2Cu306+x system. For YBa2Cu30Q, the only constraint on the hyperfine constants is the 
AF resonance result, Eq. (pj|). However, as noted above, for Y Ba2Cu^07 one has two further constraints: AB — Ac, 
and ^•^i? = 3.7 ±0.1. Moreover, as is the case for Y Ba2Cu30G,e3, neutron scattering experiments on YBa2Cu307 
suggest that one has four incommensurate and largely unresolved peaks along the zone diagonal direction whose 
positions, Qi, are given by Eq. (p^) . On taking S = 0.1, a value consistent with the experimental results of Dai et al, 



[g[ we then find, on making use of Eq. (16), that 

Aab = 0.721B. (29) 

If now we assume that the spin orbit coupling of a Cu^~^ ion in Y Ba2Cu-iO^ is little changed from that found for 
La2CuOi, 7 — 0.471, we have a third relation between the coupling constants, 

Aab -Ac^ 4.721B = 203 kOe/^lB (30) 

from which we find 

B = 43 kOe/fiB 
Aab = 31 kOe/ hb 

Ac = -172 kOe/^B (31) 



while from the AF resonance constraint, Eq. (|19[), we find for the insulator YBa2Cu30G, that B = 39.8 kOe/fiB- 

Confirmation of this choice of parameters comes by determining the slope from the linear temperature dependence 
found in a plot of ^'^Kab versus Xo{T) for Oe.es- We find 4B + Aab ^ 200 kOe/fiB, in agreement with Eq. (|3^). 
Moreover, Shimizu et al [Q find, from a similar plot for Y Ba2Cu30GAs, thai for this system, 4:B+Aab — 200 kOe/fiB- 
We adopt these values in our subsequent calculations. We note that the value of B wc obtain for Y Ba2Cu^OQ 
is some 10% larger than that found for LaCuOi, while the doping dependence of B is considerably smaller in the 
YBa2Cu306+x system than in the La2-xSrxCu04 system. Both effects may plausibly be attributed to the presence 
of chains in the YBa2Cu306+x system. The core polarization parameter, k — 0.281 we find for the YBa2Cu306+x 
system is some 10% smaller than that inferred for the La2-xSrxCu04 system. We tabulate in Table IV our results for 
the YBa2Cu306+x system at three doping levels; we estimate B = 40.6 kOe/fiB for YBa2Cu30es3 by interpolating 
between an assumed value, B = 39.8 for YBa2Cu30e.5, and that we found above for Y Ba2CuzO'j. 

IV. RECONCILING NEUTRON SCATTERING AND ''^Cu NMR MEASUREMENTS IN La2_xSrxCu04 

We now explore whether, with the revised hyperfine constants proposed above, we can reconcile the recent neutron 
scattering results of Aeppli et al for Lai^^^BrQuCuO/j, with the NMR measurements of Ohsugi et al |33 on the 



two adjacent systems, Lai.87S'ro.i3Cw04, and iai.855'ro.i5Cu04. We assume that x{Qi^) takes the MMP form, Eq. 
(^, in which case 



X Ig^^J ~ / ^ n I (r\ \2t:2n I f — 1 \2 y'^'^l 
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where XQi is given by Eg. (|To|) . There are three undetermined parameters; a, ^, and losf- We begin by deducing XQi 
and ujgp from the results of AeppU el al for x"{Qii^) at SSfsT; as may be seen in Fig. 2, a good fit to their results is 
found with xg. (SSiiT) = 350 states/eV and losf = 8.75meV. To determine a, and hence ^(SSX), we turn to the NMR 



results of Ohsugi et al; |33| on interpolating between their results for the adjacent systems, as shown in Fig. 3, we find 
63TiT = M{10-^ sK), while according to Table III, one has 



^^TiT ^M.2{ujsFla)sK/{eVf. (33) 

Equating these results, we obtain a = 23.9staies/eF and ^ = 7.6. 

A first check then on our use of Eq. (^) to fit both NMR and neutron scattering results is to compare this value of 
^ with the measurements of the intrinsic line width of each peak by Aeppli et al. We find on converting units, that 
at SSisT the linewidth parameter of Aeppli et al corresponds to a correlation length, ^ = 7.7 in the low {uj — OmeV) 
frequency limit. The agreement is quite good. 

Having determined a, we can then use our interpolated NMR results to obtain ujsf{T) for 35K < T < 300K from 
Eq. (HI). That leaves only one parameter, XQi (or to be determined over this temperature range. As a first step 
toward its determination, we use the results of Aeppli et al for x"{Qii ^) 8QK. As shown in Fig. 2, a good fit to the 
experimental data is obtained with XQii^^K) — 175 states/eV. From Eq. (^, we then get ^{80K) — 5.41. 

We next make use of the Barzy kin- Pines magnetic phase diagram. From their analysis of NMR, transport and static 
susceptibility experiments, they conclude that the La2_xSrxCu04 system will, like its YBa2Cu306+x counterpart, 
exhibit non-universal scaling behavior, perhaps best described as pseudoscaling, between two cross-over temperatures, 
T* and Tcr- In this regime, the system exhibits apparent z = 1 dynamic scaling behavior, with losf varying linearly 
with temperature and 

tosF = c'/C (34) 

where c' depends on the doping level. They propose that the upper cross-over temperature, Tcr, which marks the 
onset of pseudoscaling behavior, can be identified as the maximum in the measured value of Xo{T), and corresponds to 
a magnetic correlation length, ^ ~ 2. The lower temperature T* is determined from ^^Ti measurements as the lower 
limit of the linear variation of ujsf (or ^^TiT) with temperature. Inspection of Fig. 3 shows that for iai.86'5'?'o.i4Cu04, 
one has a comparatively weak cross-over at T* ^ 80K. Since losf{T) has already been determined, a knowledge of 
c', obtained at one temperature between T* and Tcr, enables one to fix f (T) over the entire temperature range. From 
our fit to the neutron data at 80K, we find c' = 52.9meV, and use this result to conclude that Tcr ~ 325/^, and that 

i = 0.0828 + 0.128(^) 80K<T< 325K (35) 

We can interpolate between this result for ^(T) and our result at 35K to obtain ^(T) over the region, 35K <T< 300-?^. 
The result of that interpolation, which is very nearly a continuation of the linear behavior found above 80K, is given 
in the inset of Fig. 4. 

A first check on the correctness of this procedure is to compare our "NMR" derived results at 300K, shown in 
Table II, with the neutron scattering results at this temperature. As may be seen in Fig. 2, the slope, obtained from 
the NMR results, [x"(Q,w)/w] = XQi/^SF = 1-10 ii%/{eV ■ meV) is in good agreement with experiment. A second 
check is to compare our results for ^(T) with the values deduced from half- width of the incommensurate peaks in 
Imx(q, w) observed in neutron scattering over the entire temperature domain (35 <T< 300K); that comparison is 
given in the main portion of Fig. 4. Finally, we can compare the predictions of Eqs. (jsj) and ( ^sf ) (the parameters 
being specified in Table II) with the combined frequency and temperature dependence of the half-width found by 
Aeppli et al in Fig. 5. In obtaining this figure, we calculated the theoretical inverse correlation length k from the ^{T) 
shown in Fig. 4, by matching the full width at half maximum of the incommensurate peaks of Eq.(p2|) to those of the 
experiments of Ref. |jl|. Our comparison of the calculated k{lu,T) to the experimental values is shown in Fig. 5. The 
extent of the agreement between our calculations and experiment suggests that we have succeeded in reconciling the 
^^Cu NMR results with the neutron scattering results, and it suggests as well that the neutron scattering results are 
consistent with z — \ pseudoscaling behavior for temperatures less than 300K. The latter conclusion was also reached 
by Aeppli et al from their analysis of their neutron scattering experiments. 
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V. "O RELAXATION RATES FOR Lai.85Sro.i5Cu04 



We now demonstrate that by choosing a reasonable next nearest neighbor hyperfine coupling contribution C", 
we can reconcile the incommensurate peaks in ^) with the measured NMR relaxation rates (^^TiT)^^ for 

Lai.85Sro.i5Cu04. 

In calculating (^^TicT)^^ for Lai.85Sro.i5Cu04, we simply use the previously determined parameters as inputs to 
Eqs.(||) and (^), where the next nearest neighbor Cu-oxygen hyperfine coupling C" is included in the form factor ^'^F^. 
For La2-xSrxCu04 materials, there is still not enough experimental data to determine the exact values of C^'^ for 
different field orientations; we therefore assume that these values are the same as those of the YBa2Cu306+x family. 
Following Monien et al [§ and Yoshinari et al (2|], we take Cf"^ = 33 kOe and C|| = 1-42, Cc = 1 and C± = 0.91. 
We further assume an isotropic C", with r = C /Cc = 0.25, and obtain XoiT") by modifying the results of Ref. ||] 
to reflect the new values of Aab and B presented in Sec. III. We use the a for Lai.86Sro.i4Cu04 obtained from the 
neutron scattering fits from the last section, and obtain ujsf-, and ^(T) from NMR data of Ohsugi et al. These 
numbers are almost the the same as those of Lai.86Sro.i4Cu04. The remaining parameter in Eq.(^), F, is chosen to 
get the best fit to the experimental results for {^''TicT)~^ . It is important to point out that our choice of F does not 
affect (^■^TicT)"^ and 1/^^T2g, because the Fermi Liquid contribution to these quantities is negligible compared to 
that of the anomalous antiferromagnetic spin fluctuations. 

In Fig. 6 we compare our calculated "'^^O NMR relaxation rate (^^TicT)^^, using F — 345meV, with the experimental 
data of Walstedt et al Q. The agreement is quite good. Note, however, the choice of r and F is not unique in our 
calculations; fits of the same quality can be obtained by choosing other values for r and F. The inset of Fig. 6 shows 
the substantial leakage of the anomalous spin fluctuations (the first term only in Eq(^) to (^^TicT)~^, calculated 
with the standard SMR form factor (r = 0). The Q-'^TicT)^^ thus calculated has a temperature dependence similar to 
that of the ^^Cu relaxation rates, much faster than seen experimentally. Also shown in the inset is the substantially 
smaller AF leakage calculated from the present hyperfine coupling ^^Fc with r = 0.25. 



VI. NEUTRON SCATTERING LINE WIDTHS AND ^^O SPIN-LATTICE RELAXATION RATES IN 

YBaaCusOe+x 

We turn now to the neutron scattering and NMR experiments for the YBa2Cu306+x system. As noted in the 
introduction, one apparent problem here has been that the large q-width of the antiferromagnetic peak, as observed 
in the neutron scattering experiments P-pT[, appeared to be in contradiction with size of the correlation length 

> 2) required to explain the ^^O NMR experiments. As Thelen and Pines demonstrated the half- width at 
half maximum for the antiferromagnetic peak in ^"(qi^) should have been qi/2 % 0.4/a in order to be consistent 
with the Mila-Rice-Shastry model and the oxygen relaxation data for YBa2Cu307. They found that in order to 
be consistent with experiment the leakage from the antiferromagnetic peak should account for no more than 1/3 
of the total measured oxygen rate. This upper bound from NMR is much smaller than the actual q-width of the 
antiferromagnetic peak, qi/2 — 0.7/a, observed in the neutron scattering experiments Q. Assuming the measured 
width is produced by incommensuration, we plot the antiferromagnetic "leakage" contribution (i.e., that from the 
antiferromagnetic part of Eq. (|p)) to the ^^O relaxation rate in Fig. 7, using the incommensuration 5 = O.Itt, which 
provides a fit to the neutron scattering experiments j^. Obviously, as in the La2-xSrxCu04 material, the temperature 
dependence of the measured NMR relaxation rate is remarkably different, and the amplitude of the "leakage" term 
is too large. This problem can be avoided by introducing C", as we have done on the La2_xSrxCu04 system. In fact, 
the much smaller degree of presumed incommensurability in the YBa2Cu306+x system than that measured directly 
for the La2-xSrxCu04 system makes it almost evident that any problem produced by AF leakage can be reconciled 
by the same method as used above. We show, in Fig. 7, that the AF leakage contribution for r — C /Cc — 0.25 indeed 
becomes negligible. If we assume the same ratio of the AF part to the total rate as Thelcn and Pines did, we 
obtain a constraint on C". We note that the oxygen form factors Eq.(||) are quadratic in (5q = (q — Q) in the vicinity 
of the antiferromagnetic wave vector Q: 

i7F„ = (5g.f-L^J— [(2rc + l)Ca,-4r„J^=rK<5<z.f (36) 

ai—a' .a" 

As a result, the antiferromagnetic contribution to the oxygen relaxation rate (which we keep as a constant when we 
change the form factor) is 



9 



1 



«77(Qi-Q)', (37) 

AF 



and a change of the oxygen form factor, which alters 77, produces a new constraint on the acceptable width (or 
incommensurability) of the neutron scattering peak. Since Thelen and Pines pO|] used the isotropic form of the 
Mila-Rice-Shastry Hamiltonian, with Ciso ~ Cf'^^ we easily obtain from Eq.(|3^): 

gi/2a < 0.4/a , ^^^'^^ (38) 

1 [(2r, + l)C., -4r„J^ 



where we have neglected possible slow (logarithmic) dependence. In particular, with — 0.25, Eq.(38) gives the 
upper limit: qi/2 ^ 0.7/a. This crude estimate shows that indeed, our hyperfine Hamiltonian is consistent with 
both NMR and neutron scattering experiments. However, the antiferromagnetic leakage contribution to the oxygen 
relaxation rate in YBa2Cu306+x can become important, and should therefore be calculated numerically, since the 
spin- spin correlation length is very short. 

For our numerical calculation of the antiferromagnetic peak contribution to the ^^O relaxation rates we assume, 
as indicated in the Introduction, that the neutron scattering data of Tranquada et al |pl| and Dai et al ^ can be 
interpreted as indicating that the magnetic response function x(Qj ^) possesses four incommensurate peaks located at 
Qi = (7r±5, 7r±(5), and take 5 ~ O.Itt, an incommensuration consistent with the measured experimental widths. We also 
assume that the temperature-dependent spectral weight for these incommensurate peaks, as in case of La2_xSrxCu04, 
comes from the temperature dependence of the correlation length, and adopt the MMP form Eq.(^ for each of the 
four peaks. It should be emphasized, however, that accord between the inelastic neutron scattering and the oxygen 
NMR can be reached for any bell-shaped curve for x"(q, f^) which has the characteristic width measured in the 
neutron scattering experiments, and a sufficiently abrupt fall-off at large (q — Q). In Fig. 7, we show our calculated 
antiferromagnetic leakage to the oxygen relaxation ^"^Wic/T = 1.5(^^TicT)~^, for the case of both r = and r = 0.25; 
again, we see that the new form factor with r = 0.25 greatly reduces the AF leakage. Also shown in Fig. 7 is our 
calculated ^^Wic/T plotted against the experimental data of Martindale et al. In obtaining our theoretical result, 
we have used as an input to Eq.(|9|), XoiT) deduced from the Knight shift Kc{T) data on the same sample, provided 
by Martindale et al and used the ^ and a from Ref. Again, we take C°''^ = 33 kOe/^B for YBa2Cu306+x 
system. By assuming ~ 0.25, we obtain a good fit to the experimental data with F ^ 308mey. These parameters 
are listed in Table IV. 

Another problem with the one-component Shastry-Mila-Rice picture has been pointed out recently by Martindale 
1^, who measured planar ^^O relaxation rates for different magnetic field directions. They have found that the 
temperature dependences of the relaxation rates for magnetic fields parallel and perpendicular to the Cu-0 bond axis 
directions were different, in contradiction with the predictions based on the MSR hyperfine Hamiltonian for which 
the oxygen form factor is given by Eq.(^), without C": 

"^^.=2 ^os^'-f^l- (39) 

ai — a' .a" 

From Eq.(p9|) it follows that the ratios of the oxygen relaxation rates for different magnetic field orientations should 
be temperature-independent, and determined only by the Hyperfine C-couplings: 

Experimentally, as shown by Martindale et al, Q these ratios turn out to be mildly temperature-dependent, although 
numerically close to the values of Eq. (|^ . 

This apparent contradiction can, in fact, be turned into a proof of the validity of the modified one-component 
model Eq.(|l|). It can easily be seen that for our oxygen form factors, Eq.(^, the ^^O relaxation rates for different field 
directions do not have the same q-dependence for the whole Brillouin zone. As a result, ratios such as Eq.(^) should 
indeed become temperature-dependent. Since we do not know the precise values of the couplings once we go beyond 
the nearest-neighbor Mila-Rice-Shastry approximation, we use here the expressions for the oxygen form-factors in the 
most general form. To derive the form of the temperature dependence, we separate the antiferromagnetic and the 
Fermi hquid or short wave-length (xo) contributions to (1/^''Tiq), according to Eq.(^): 



1 / 1 \ / I 



(41) 

AF 
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Here the short wave-length part {l/^'^TiT)^^ is proportional to the bulk magnetic susceptibility Xo{T), while the 
antiferromagnetic part follows the copper relaxation rate: 



1 



_J_\ ^ ^^K(Qi) f 1 



Xo 

17; 



(42) 



As we have demonstrated above, the temperature-dependence of the antiferromagnetic leakage term is very different 
from what is observed in experiment. Since the empirical modified Korringa law {l/^'^TiTxo{T)) = const is rather 
well satisfied for these materials, the short wave-length part should be dominant. Therefore, we can write for the 
different ^^O relaxation rate ratios: 



''J 



S., Sj J \^^Ti,Txo{T) 



(43) 



where Pj = Fa - (Qi)/ -F'c(Qi), while the Sj are coefficients determined by integrating the product of the short-range 
part of the magnetic susceptibility with the oxygen form factor. If the short-range part of x"(q, w) is only mildly 
q-dependent, Sj is determined primarily by the momentum average of ^'^ Faj, 

S,^^ J ''FaM)d^<l (44) 

In this case the temperature-independent part of the ratio of the oxygen relaxation rates is determined again only by 
the ratio of the form factors: 

^ ^ ^"fa.(r = 0) 
Sj i7F„^(r = 0) 

If a realistic band-structure q-dependence of x"(q, a;) is taken into account, this ratio will have a somewhat different 
value. We demonstrate, in Fig. 8, that expression Eq.(^) indeed provides a consistent explanation of the temperature- 
dependent term for the oxygen relaxation rates in YBa2Cu307; on using Eq.(^3|) to fit the observed anisotropy ratio 
of we find 



S\ 



s± 

Sn /Pn P 



T^- = 0.5, 



iS* I V S*!! S_i 



0.06 {sK)i?g/eV (46) 



These values of Si_j and Pij impose certain constraints on the parametric space of the hyperfine couplings. However, 
there are not enough of these constraints to enable us to deduce unambiguously the values of the hyperfine couplings, 
so that specific quantum chemical calculations are needed to determine the hyperfine coupling constants for these 
materials. However, as we have shown, the temperature dependence of the rates can be accounted for by assuming a 
finite incommensurability for the antiferromagnetic peak. 

Using our formalism, and the constants 5*^ and Pj for YBa2Cu307, we can predict the oxygen relaxation rates for 
YBa2Cu306.63- It is easy to see that if the hyperfine C-couplings do not depend significantly on doping, the product, 
SjT, for YBa2Cu3 06.63 is the same as for YBa2Cu307. Pj^ however, can be somewhat different, corresponding to a 
different amount of incommensuration for YBa2Cu3 06.63- Since the oxygen form factor is quadratic in the vicinity of 
(7r/a,7r/a), we can write: 

p _ p "YBa2Cu306.63 I Anx 

-f^YBa2Cu306.63 - -f^YBasCuaOT-rt V^'j 

"YBaaCuaOr 

However, the degree of incommensurability in YBa2Cu3 6.63 [0 is roughly the same as in YBa2Cu307, both have 
5 ~ O.Itt, so that Pj will remain unchanged from the YBa2Cu307 values. This makes it possible to predict the 
behavior of the oxygen relaxation rates ratios in YBa2Cu3 06.63j once the parameter T in Eq.(^) is determined from 
experiment. In Fig. 9, we fit the ^''O relaxation rates (^^TicT)"^ of Takigawa et al to determine Sc- Again, we use the 
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XOi 01 and ^ given in Ref. ||]. It is seen in the main portion of Fig. 9 that the fit is very satisfactory, from this fit we 
obtain F = 22QineV for YBa2Cu3 06.63: if we assume = 0.25. These parameters are also listed in Table IV. From 
Eq.(Q, we have SiT being the same in YBa2Cu307 and YBa2Cu306.63, because their form factors do not change. 
Therefore, we get for YBa2Cu3 06.63, 

^^f^l-E^] 'gc(YBa2Cu307) 

S± \Sil S'±/YBa2Cu306.63 *^^/YBa2Cu30. S'c(YBa2Cu3 06.63) 

^ Q QQ ^ rYBa2Cu306.a3 ^ q q^^ {sK)tll/eV (48) 
i YBaaCuaO- 

We show our calculated relaxation rate ratio for YBa2Cu3 06.63 in the inset of Fig. 9. 

VII. DISCUSSION AND CONCLUSIONS 

We have seen that by modifying the SMR hyperfine Hamiltonian we can use the MMP one-component spin-spin 
response function to reconcile the results of a number of neutron scattering and NMR experiments on the cuprate 
superconductors. With the aid of the scaling arguments of Barzykin and Pines, ^ we are able to obtain a quantitative 
fit to both the NMR and the neutron scattering data for Lai.86Sro.i4Cu04. We find that for the YBa2Cu306+x system, 
we can reconcile the q- width of the antiferromagnetic peak seen in neutron scattering experiments with the substantial 
temperature dependent AF correlation required to explain the NMR experiments on YBa2Cu307 and YBa2Cu3 6.63- 
Moreover, in the recent results of Martindale et al |^ on the anomalous temperature dependence of the anisotropy of 
the ^^O relaxation rates, the small amount of the AF leakage is shown not only to be explicable using our modified 
one-component description; but to provide a direct proof for the one-component picture. Our ability to reconcile so 
many different experiments leads us to conclude that a transferred hyperfine coupling between next nearest neighbor 
Cu^"*" spins and ^^O nuclei spin plays a significant role, and that the transferred hyperfine coupling B, changes as one 
goes from the La2-xSrxCu04 to the YBa2Cu306+x system, and is moreover comparatively sensitive to hole doping 
in the former system. It will be interesting to see whether the presence of these terms can be justified microscopically 
through detailed quantum chemical calculations in these systems. 

Our results have a number of interesting implications for NAFL (nearly antiferromagnetic Fermi liquid theory) 
calculations of other properties of the superconducting cuprates. For example. Pines and Monthoux |]35[] have shown 
that incommensuration acts to lower the superconducting transition temperature, T^, it is tempting therefore to 
attribute much of the substantially difference in Tc found for the La2-xSrxCu04 and YBa2Cu306+x systems to the 
much greater degree of incommensuration found in the former materials. In their calculation of planar resistivities, 
Stojkovic and Pines find that pab depends sensitively on the size and distribution of "hot spots" ( regions of 
the Fermi surface connected by Qi), and thus is markedly changed by incommensuration. To cite a third example, 
in NAFL theory, the location in momentum space of the peak in the spin fluctuation spectrum depends on the 
interplay of the peaks in the irreducible particle-hole susceptibility, x(q, 0), produced by band structure and the 
momentum-dependence of the restoring force, Jq, which acts to shift those peaks according to Ref. ^T\ , 

1 - Jqx(q,o) 

Since the peaks in xiq.^ 0) move away from (tt/o, tt/o) as one moves away from half-filling, less peaking in Jq is required 
to produce four incommensurate peaks than was needed by Monthoux and Pines ||3^ to keep the peak at {tt / a^ir / a) 
in the presence of substantial hole doping. 

Further NMR and neutron experiments on the YBa2Cu306+x and La2-xSrxCu04 systems can also help verify 
the correctness of our proposed new hyperfine Hamiltonian and our assignment of incommensurate peaks in the 
YBa2Cu306+x system. For example, our results, Eq.(|36|) and Eq.(^3|), lead us to predict substantial temperature 
dependence in the anisotropy of 1/^^Tiq in the La2_xSrxCu04 system for magnetic fields parallel and perpendicular 
to the Cu-0 bond axis, and it will be instructive to see whether this can be measured. It is, moreover, to be hoped 
that improvements both in neutron scattering facilities and the availability of large single crystals will make possible a 
direct experimental check on our assignment of incommensuration in the YBa2Cu306+x system. Resolution of those 
peaks, together with a direct measurement of their intensities would also enable one to carry out a detailed comparison 
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of NMR and neutron scattering experiments on YBa2Cu306.63 analogous to that presented here for Lai.86Sro.i4Cu04 
system. 
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TABLE I. Spin-lattice anisotropy and incommensuration in the La2-xSrxCu04 system 



System 


S 


J^Jexpt 


Ref. 




-K_Bq.(14) 


LaC'u204 
Lai.gSro.iCu204 
Lai,s5Sro.i5Cu204 



.175 
.263 


3.9±0.3 
3.5±? 
3.0 ±0.20. 




3( 
~ 

31 


1 
1 
1 


3.7 
4.11 
4.78 


3.9 
3.2 
3.2 



TABLE II. Fits to neutron scattering experiments of Lai.86Sro.i4Cu04 





T=35K 


T=80K 


T=325K 


''''TiT(xlO-'sK) 
<^SF(meV) 
a(states/eV) 
XQ(Ml/eV) 

e (a) 
c' 


34 
8.75 
23.9 
350 

7.6 


38 
9.78 
23.9 
175 
5.41 
52.9 


103 
24.7 
23.9 
26.2 
2.14 


l/T2G(msec-^) 


63 


45 


18 



TABLE III. Parameters for La2_xSrxCu04 





La2Cu04 


Lai.9()Sro.ioCu04 


Lai.86Sro.i4Cu04 


Lai.85Sro.i5Cu04 


Ac(kOe/^is) 
Aab(kOe//iB) 
B(kOe/^fl) 
Cc(kOe/^ifl) 

63 R 

J^cxp 
63r 

4B-Aab(kOe/^is) 
4B+Aab(kOe//is) 

e^'Ke/^^^Kab 

63 r-|-i^ r-|-i 

1/T2G 
63TiT/T2G 

^^TiT/Tic 


-185 

18 
36.1 

33 
3.9±0.3 

3.9 

127 
162.4 
-25% 

l3S{sR'/eV'^)ujsF/a 

298(eV/s)aC 
i.l2xW^{K/eV)LosF^ 
1.23x10'' {K/s)aujsF^^ 


-185 
18 
46 

33 
3.5±? 
3.2 
166 
202 
-0.5% 
95.2{sK/eV^)ujsF/a 

347(eV/s)QC 
3.30xlO^{K/eV)ujsF^ 
1.15x10'^ {K/s)aujsF^^ 


-185 
18 
50 
33 

3.2 

182 
218 

7% 

93.5{sK/eV'^)LOsF/a 

350(eV/s)aC 
3.27xlO'^{K/eV)LosFi 
l.UxlO''{K/s)aujsF^^ 


-185 
18 
51 
33 
3.0±0.2 
3.2 
186 
222 
8.6% 

9'i.2{sK/eV^)uJsF/a 

348(eV/s)aC 
3.28xlO'^{K/eV)LosF^ 
l.UxlO\K/s)aujsF£,^ 


r (meV) 
r 
6 




0.175 


0.245 


345 
0.25 
0.263 
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TABLE IV. Parameters for YBa2Cu306+x 





YBa2Cu306 


YBa2Cu306.63 


YBa2Cu307 


Ae(kOe/^is) 

Aab(kOo/^B) 

B(kOe/Ms) 
Cc(kOe/Ms) 

63 R 
J^exp 


-172 

31 
39.8 
33 

3.8 


-172 
31 

40.6 
33 

4.0 


-172 
31 
43 
33 
3.7±0.1 
3.7 


4B-Aab(kOe/MB) 
4B-t-Aab(kOe//XB) 

63r|-i^ rj-i 

I/T2G 

ill/ i2G 

"^^TiT/TIg 


128.5 
190 
-7% 

135(s7f/el/^)tJsF/a 
301(eV/sX 

4i:.UDXlU l^i\/eK jiOsFC, 

1.22x10'^ {K/s)auJsFi^ 


131.4 
193 
-5% 

U5isK/eV^)uJsF/a 
293(eV/sX 

1.25x10'' {K/s)auJsFi'^ 


141 
203 


126{sK/eV'^)LOSF/a 

310(eV/s)aC 
o.yxiu \^j\/cv )iOsF<, 
1.21x10' {K/s)auJsFi^ 


a 




8.34 


15.36 


r (meV) 




226 


308 


r 




0.25 


0.25 


5 




0.1 


0.1 



15 




16 




FIG. 2. The frequency dependence of x"iQi,^) at three temperatures. The experimental points are the results obtained by 
Aeppli et al the solid curves are the fits obtained using a mean field description of x"(q, t^) shown in Eq.(^^ and parameters 
compatible with NMR results. 
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FIG. 3. The interpolated TicT for Lai.86Sro.i4Cu04 is shown together with the measured values of TicT for 
Lai.g7Sro.i3Cu04 and Lai.85Sro.i5Cu04 of Ohsugi et al. js^ Shown on the right hand side is the scale for ujsf{T) oc ^^TicT 
for Lai,86Sro.i4Cu04 inferred from the fit to the neutron scattering experiments. 
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FIG. 4. A comparison of the NMR-deduced values of ^(T) (solid line) for Lai.86Sro.i4Cu04 with those obtained (diamonds) 
from the neutron scattering experiments of Ref. [Q]. The experimental points (diamonds) are derived by first fitting the 
half- width of the neutron scattering peak at Qi for low energy transfer {uj = 2.5meV), and then extrapolating to w = 0, 
following the formula K^(tj,T) = k^(T) + a~^Lu^/E^ given in Ref. [Q]. The inset shows the interpolation procedure used to 
obtain 1/C(T) between 35K and 80K: the point at 35K (star) is obtained from the MMP fit to the neutron scattering data at 
35K, while the points above 80K are deduced from the scaling analysis of Eq.(^) (circles); the solid line shows the extrapolation 
between 35K and 80K. 



19 



0.20 




0.00 I — ' — ' — ' — 

1 10 100 1000 

T(K) 



FIG. 5. A comparison of the NMR-deduced values of the frequency dependence inverse correlations length, k, at tj=3.5, 
6.1, and 15meV (lines), with the experimental results of Aeppli et al |l| (symbols). It is seen that the consistency is quite 
good at low frequencies, while at high frequencies, the NMR-deduced q-width is smaller than those seen in neutron scattering 
experiments. 
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FIG. 6. Our calculated spin-lattice relaxation rates (^^TicT)^^ for ^^O (solid line) compared to the experimental data of 
Walstedt et al |^ (circles) for Lai.85Sro.i5Cu04; we use r=0.25 for the calculation. The inset shows the contribution of the 
AF spin fluctuations (first term only in Eq.(^)) to (^^TicT)"^, calculated with the present form factor (r=0.25) and with the 
standard Shastry-Mila-Rice form factor (r=0). 
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FIG. 7. The spin-lattice relaxation rate ^^Wic/T calculated by assuming r=0.25 (solid line), plotted against the ex- 
perimental data of Martindale et al ^ (circles) for YBa2Cu3 06.96- Also shown is the contribution from AF leakage to the 
relaxation rates ^"^Wic/T calculated using our oxygen form factor with r=0.25 and the standard Shastry-Mila-Rice form factor 
(r=0). 
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FIG. 8. The temperature dependence of the oxygen relaxation rate ratios in YBa2Cu3 06.96 measured by Martindale et al 
|, and fits using our theoretical expression Eq([43|). 
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FIG. 9. The calculated spin-lattice relaxation rates C^TicT)~^ for YBa2Cu3 06.63, compared with the data of Takigawa 
et al. [^] The inset shows the predicted anisotropy ratio as described in the text. The predicted anisotropy lies slightly above 
the calculated curve for YBa2Cu3 06.96 in Fig. 8, yet both are within the experimental error bars of the YBa2Cu3 06.96 data of 
Martindale et al. Q 
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